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ABSTRACT 

We have detected the optical counterpart to the nearby isolated neutron star RX J1605. 3+3249 using 
observations from the Space Telescope Imaging Spectrometer aboard the Hubble Space Telescope. The 
counterpart, with msocCD = 26.84 ± 0.07 mag and very blue colors, lies close to the ROSAT HRI error 
circle and within the Chandra error circle. The spectrum is consistent with a Rayleigh- Jeans tail whose 
emission is a factor of w 14 above the extrapolation of the X-ray blackbody, and the source has an 
unabsorbed X-ray-to-optical flux ratio of log(/x//opt) = 4.4, similar to that of other isolated neutron 
stars. This confirms the classification of RX JI605. 3+3249 as a neutron star. 

Subject headings: pulsars: individual (RX J1605. 3+3249) — stars: neutron — X-rays: stars 



1. INTRODUCTION 

Thanks to ROSAT, a half dozen nearby neutron stars 
that emit no detectable radio emission have been identified 
(see Treves et al. 2000 for a comprehensive review) . These 
objects have been eagerly studied using facilities such as 
the XMM-Newton Observatory, the Chandra X-ray Ob- 
servatory, the Hubble Space Telescope (HST), Keck and 
the Very Large Telescope with the hope of determining 
the physical parameters, in particular radius and temper- 
ature, and to compare these to models of neutron stars. 
Independently, by their sheer proximity these objects play 
a pivotal role in assessing the neutron star demographics in 
the Galaxy. These two considerations highlight the virtue 
of detailed studies of nearby neutron stars. 

RX J1605.3+3249 was identified in the ROSAT All-Sky 
Survey by Motch et al. (1999). The X-ray spectrum is 
well fitted by a blackbody with kT ~ 90 eV and low in- 
terstellar column density, Njj ~ 10 20 cm~ 2 , and is quite 
similar to those of well studied nearby neutron stars such 
as RX JI856.5-3754 and RX J0720.4-3125. Motch et al. 
(1999) obtained deep (B ~ 27 mag and R ~ 26 mag) im- 
ages from the Keck telescope. Only one object (star C) 
was found within the 2" High Resolution Imager (HRI) 
circle (we believe that the uncertainty of the HRI position 
was underestimated; see § 2.2). Optical spectroscopic ob- 
servations carried out at the Canada-France-Hawaii Tele- 
scope (CFHT) showed that star C was a distant late-type 
M dwarf. 

However, the soft spectrum and the stable X-ray 
emission are better accounted for by a model in which 
RX J1605. 3+3249 is an isolated neutron star. If so, the 
optical counterpart would be below (or perhaps just at) 
the limit of the Keck observations. As a part of our inves- 
tigation of nearby neutron stars with HST we undertook 
deep observations of this field. In this Letter we report 
the discovery of a faint blue optical star which we iden- 
tify with the optical counterpart of RX J1605. 3+3249. 
Our discovery confirms that RX J1605. 3+3249 is a nearby 



neutron star. 

2. OBSERVATIONS & DATA REDUCTION 

2.1. Hubble Space Telescope Observations 

We observed RX J1605. 3+3249 with the Space Tele- 
scope Imaging Spectrograph (STIS) aboard HST in two 
modes: unfiltered CCD (50CCD aperture) and a longpass 
filter that transmitted photons longward of « 5500 A; see 
Table 1. For each mode, the individual images were driz- 
zled (Fruchter & Hook 2002) onto a single image with a 
pixel scale of 0.5. Thus final images had 0'.'0254 pixels. 

2.2. Keck Observations 

For astrometric purposes we obtained imaging data from 
the Low- Resolution Imaging Spectrometer (LRIS; Oke 
et al. 1995) on the Keck II telescope; see Table 1. These 
images were reduced in a standard manner using tasks in 
MIDAS: bias subtraction, flat-fielding, and stacking of the 
exposures. We show the LRIS image in Figure 1. 

After correcting the stellar positions measured in a 30-s 
LRIS image for geometric distortions 1 , we fit for plate- 
scale, zero-point, and rotation using 20 unsaturated stars 
from the latest version of the Guide Star Catalog (GSC- 
2.2) 2 , obtaining rms residuals of 0'.'09 in each coordi- 
nate (in what follows, astrometric uncertainties refer to 
rms values in each coordinate unless otherwise specified). 
We then determined the astrometric solution (plate-scale, 
zero-point, and rotation) of the full 30-minute LRIS image 
by using 24 stars common to that and the 30-s image, get- 
ting rms residuals of 0'.'02. Finally, we used 21 stars on the 
deep LRIS image to determine the plate-scale, zero-point, 
and rotation of the drizzled STIS 50CCD image (shown 
in Fig. 2), getting rms residuals of 0'.'04. The final uncer- 
tainty to which our STIS coordinates are on the ICRS is 
dominated by the 0'.'3 uncertainty of the GSC-2.2 3 . 

1 See http:/ /alamoana. keck. hawaii.edu/inst/lris/coordinates. html. 

2 See http://www-gsss.stsci.edu/support/data_access.htm. 

3 See http:/ /www-gsss.stsci.edu/gsc/gsc2/calibrations/astrometry/astrometry.l 
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Table 1 

Summary of Optical Observations 



Telescope 


Instrument 


Date 


Exposure 


Band 






(UT) 


(s) 




Keck II 


LRIS 


1998-Aug-24 


1800 


R 


HST 


STIS 


2001-Jul-21 


2700 


50CCD 


HST 


STIS 


2001-Jul-21 


5360 


F28X50LP 



LRIS R-band 




Fig. 1. — LRIS R-band image of the field around 
RX J1605.3+3249. The field-of-view of the STIS 50CCD obser- 
vation is indicated by the box. The 15"-radius (90% confidence) 
PSPC and l'.'5-radius (90% confidence) updated HRI error circles 
are also indicated. Sources A, B, C, and D from Motch ct al. (1999) 
arc labeled. 



We inspected all point sources on the STIS images that 
were within the PSPC error circle and found a very blue 
source. As can be seen from Figure 3 this is the bluest 
source in the PSPC error circle. The source, hereafter "X" , 
is located at (J2000) a = 16 h 05 m 18?52 S = +32°49'18"0, 
with uncertainties of about 0'.'3, and lies 2'.'5 from the HRI 
position from Motch et al. (1999), outside the nominal 90% 
error circle. 

However, we have found a problem in comparing posi- 
tions referenced to the USNO-A2.0 (Monet 1998), like the 
HRI position, and those referenced to the GSC-2.2/ICRS 
(our optical data). In this field, there appears to be a sys- 
tematic shift between the GSC-2.2 and the USNO-A2.0 of 
(<5usno — <5gsc) = 0'.'65 (the shift in Right Ascension is 
a negligible Of! 02). If we correct the HRI position to the 
GSC-2.2 reference frame, we find the new position to be 
(J2000) a = 16 h 05 m 18!66 S = +32 o 49'19'.'0. With this po- 
sition, the error circle appears to be located properly with 
respect to star C, comparing with Fig. 4 of Motch et al. 
(1999). 

Even with the updated position, X is slightly outside the 
90% HRI error circle 4 . We have two possible explanations 
for this. First, RX J1605. 3+3249 may have non-negligible 

4 X is within the error circle of the position obtained from a prelim- 
inary analysis of XMM data; van Kerkwijk et al., in preparation. 



proper motion, such as that seen for RX J 1856. 5— 3754 
(Walter 2001). In this case the offset between the HRI po- 
sition (epoch 1998.3) and the STIS image (2001.6) could 
be real. However, the LRIS data are not of sufficient qual- 
ity to detect X with any confidence, so we will have to 
wait for additional data. Second, we note that the quoted 
uncertainty of the HRI position may be underestimated: 
Motch et al. (1999) used 6 reference sources for the bore- 
sight corrections and claim an uncertainty of Of! 64 with 
no contribution from systematic effects. In comparison, 
Hasinger et al. (1998) use 32 sources and get typical HRI 
uncertainties of 1.0" that include a 0"5 systematic error to 
achieve good X-ray-to-optical matches. 

In either case, the blue color of source X is similar 
to those of the counterparts of other isolated neutron 
stars (e.g. RX J1856. 5-3754 and RX J0720.4-3125; Wal- 
ter & Matthews 1997; van Kerkwijk & Kulkarni 2001b; 
Motch & Haberl 1998; Kulkarni & van Kerkwijk 1998). 
Thus we consider it likely that X is the counterpart of 
RX J1605. 3+3249. 

2.3. X-ray 

While our identification based on color and position is 
plausible, the uncertainty in the ROSAT HRI position pre- 
vents us from being sure about the association (pulsations 
or a common proper motion would be definitive measure- 
ments). Fortunately, the availability of archival data from 
the Chandra X-ray Observatory offered us to opportunity 
to decrease the chance coincidence probability by a fac- 
tor of 10. The Chandra observation (ObsID 2791) had a 
duration of 20-ks, and RX J1605. 3+3249 was at the aim- 
point of the ACIS-I CCD array. Using standard processing 
steps 5 we corrected for a systematic astrometric error of 
Aa = — 0"35 and AS = —Of! 10. As a cross-check, we com- 
pared the positions of other X-ray sources to GSC-2.2 stars 
(which, since we used the GSC-2.2 as the reference for the 
optical astrometry, ensures that they are on the same sys- 
tem as our HST data) and found that the coordinates 
match to better than 0.5 arcsecond. 

We then measured the centroid of RX J1605.3+3249 
(with a count-rate of w 0.15 s -1 , RX J1605.3+3249 is 
somewhat affected by photon pileup, but this should not 
affect the centroid) to be (J2000) a = 16 h 05 m 18!50, 
S = +32°49'17"4. We estimate a final 90% confidence 
radius of the X-ray position with respect to the STIS im- 
age of « l'.'O. As can be seen in Fig. 2, source X is well 
within this radius, lending credence to our identification. 

5 http://asc.harvard.edu/cal/ASPECT/fix_offset/fix_offset.cgi. 
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Fig. 2.— STIS/50CCD image of RX J1605. 3+3249. The image 
has North up and East to the left. A 5" scale bar is in the upper 
right. The updated 90% confidence error circles from ROSAT /HRI 
(1''5 radius) and Chandra (l'.'O radius) are plotted. The likely coun- 
terpart to RX J1605. 3+3249 is indicated by the X, and sources B 
and C from Motch et al. (1999) are labeled (source B is not just 
extended to the South-East, as noted by Motch et al. 1999, but is in 
fact composed of separate sources). The PSPC error circle (Fig. 1) 
is larger than this image. 



3. ANALYSIS & DISCUSSION 

We rely on the spectral fits to the ROSAT PSPC data 
presented in Motch et al. (1999): kT = 92 eV, N H = 1.1 x 
10 20 cm~ 2 , and = 3.3d 30 o, where d — 300^300 pc is the 
distance and the normalization assumes 0.9 counts s _1 in 
the PSPC. The absorption column density implies an ex- 
tinction of A v = 0.06 mag, (Predehl & Schmitt 1995) 
- while this value is uncertain, it is low enough to not 
make a large difference. The total Galactic hydrogen col- 
umn density is 2.4 x 10 20 cm~ 2 (determined by COLDEN 6 ; 
Dickey & Lockman 1990) so that the maximum extinc- 
tion is 0.13 mag. This agrees with the extinction esti- 
mated from infrared dust emission (Schlegel, Finkbeiner, 
& Davis 1998), and confirms that the total extinction to 
RX J1605. 3+3249 is low. 

Within a O'.'25-radius aperture, we measure magnitudes 
of m 50 ccD = 27.03 ± 0.07 mag and to F 28X50LP = 28.16 ± 
0.18 mag for source X in the STMAG system. To cor- 
rect the photometry to a nominal infinite aperture, we 
follow Kaplan, Kulkarni, & van Kerkwijk (2002a) and Ka- 
plan et al. (2003). As X is bluer than all of the stars 
in the image, we used the bluest of the available aper- 
ture corrections: 0.183 mag at 0'.'25 radius for 50CCD and 
0.214 mag at 0"25 radius for F28X50LP (T. Brown 2002, 
private communication). This then gives corrected mag- 
nitudes of m5occD = 26.84 ± 0.07 mag and mF28X50LP = 
27.95 ± 0.18 mag, where we have incorporated a 0.02 mag 
uncertainty from the aperture corrections. 

Since the STIS bandpasses are so wide we must use the 
shapes of the bandpasses to convert the measured mag- 
nitudes into fluxes. Following van Kerkwijk & Kulkarni 
(2001b), we find effective wavelengths (A) of 5148 A and 
7137 A and effective extinctions (A x /A v ) of 1.56 and 0.79 

6 See http://asc.harvard.edu/toolkit/coldcn.jsp. 



for 50CCD and F28X50LP, respectively. With these wave- 
lengths, we can now apply the standard STMAG conver- 
sion of F\((X)) = icr^ 21 - 1 )/ 2 - 5 ergs s" 1 cuT 2 A" 1 . 

We plot the spectral energy distribution of source X in 
Figure 4, assuming a blackbody spectrum in the X-ray 
regime 7 . As one can see, the optical photometry appear 
to follow a power-law with spectral index a « 2 {F u ccv a ), 
appropriate for a Rayleigh-Jeans tail and similar to that 
of other isolated neutron stars (van Kerkwijk & Kulka- 
rni 2001b; Kaplan et al. 2003). We also plot in Figure 4 
the extrapolation of the best-fit ROSAT PSPC blackbody. 
The Raylcigh- Jeans fit to the STIS data has a normaliza- 
tion that is a factor of 14 ± 2 above the extrapolation of 
the blackbody. Other power-law indices are possible, as 
found for RX J0720.4-3125 (Kaplan et al. 2003), but we 
do not believe that the current data warrant a full fit. 
The unabsorbed bolomctric X-ray-to-optical flux ratio is 
log(/x// pt) = 4.4 (assuming that the X-ray spectrum is 
well described by a blackbody). 

The blueness of source X and its X-ray-to-optical 
flux ratio, taken together with the X-ray spectrum 
of RX J1605. 3+3249, virtually guarantee that X is 
the optical counterpart of RX J1605. 3+3249 and that 
RX J1605. 3+3249 is an isolated neutron star. It then joins 
3 other sources (RX J1856.5-3754, RX J0720.4-3125, and 
possibly RX J1308.6+2127) that have soft X-ray black- 
bodies, blue optical counterparts, and no other emission; 
see Table 2. RX J1605. 3+3249 stands out from the other 
sources in Table 2 by virtue of its relatively large opti- 
cal excess: the optical flux is a factor of 14 above the 
extrapolated X-ray flux, where for the other sources the 
ratio is closer to 6 (this is despite the fact that the X- 
ray-to-optical flux ratio for RX J1605. 3+3249 is within 

7 This blackbody is likely a simplification of a more realistic at- 
mosphere model, but as yet such models have been unsuccessful in 
fitting sources like RX J1856.5-3754 and RX J0720.4-3125. 
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m 50CCD- m F 2 8X50LP< S ™ AG > 

Fig. 3. — Color-magnitude diagram of point-like source from 
the STIS data. 50CCD magnitude is plotted against 50CCD - 
F28X50LP color. Source X is indicated and plotted as a square. 
Source C is also indicated. The sources have been corrected 
to "infinite" aperture using the correction most appropriate for 
RX J1605.3+3249 (§ 3). A reddening vector is plotted for A v = 
0.1 mag. 
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Fig. 4. — Spectral energy distribution for source X corrected for 
absorption with Ay = 0.06 mag. The STIS data arc plotted as 
points. The extrapolation of the ROSAT blackbody (Motch et al. 
1999) is the solid line (labeled "ROSAT BB"), and a Ray leigh- Jeans 
fit to the STIS data is the dashed line. The horizontal error-bars 
show the bandpasses of the filters. 



the range of the other sources). This can also be seen 
from the ROSAT data in Schwope et al. (1999), where 
RX J1605. 3+3249 has a factor of 2 smaller count rate 
than RX J0720.4— 3125 despite being hotter and having a 
comparable optical magnitude. There are two possibilities 
to explain the large excess of RX J1605. 3+3249: either its 
X-ray emission is suppressed relative to the optical, or the 
optical emission is enhanced. 

The hrst scenario implies that the blackbody radius 
(3.3<i3oo km, where the distance of <~ 300 pc is predicted 
based on the optical flux; see Kaplan, van Kerkwijk, & 
Anderson 2002b) is significantly smaller than those of the 
other sources (typically 6 km) . If the blackbody radius can 
be interpreted as the radius of a hot polar cap, perhaps 
RX J1605. 3+3249 has a different magnetic field configura- 
tion leading to a smaller cap size, or RX J1605. 3+3249 is 
in an orientation where only half of the cap is visible (al- 
though this is difficult when relativistic beaming is taken 
into account; Psaltis, Ozel, & DeDeo 2000). 

The second scenario could occur if there were a sig- 
nificant contribution to the optical emission from non- 
thermal emission. Non-thermal emission could arise if 
there were a substantial spin-down luminosity, E (such 
as that seen for PSR B0656+14; Koptsevich et al. 2001). 
In this case, the high optical excess could indicate a large 
E for RX J1605. 3+3249. While we cannot constrain the 
non-thermal emission from the current photometry, it may 
be difficult to reproduce the thermal-like spectrum ob- 
served in Figure 4 and to invoke significant non-thermal 
flux. One could appeal to light-element atmospheres for 
RX J1605. 3+3249, as they have high optical excesses over 
blackbodies (Pavlov et al. 1996). However, these excesses 
are far too high (factors of 50-100, instead of 14) to fit our 
photometry, allowing us to reject such models. 

Future observations, such as higher-precision X-ray 
spectroscopy from Chandra and XMM (being analyzed), 
additional optical photometry, and improved X-ray timing 
(in order to determine P and eventually E and the mag- 



netic field B), should help to settle these issues. We are 
also searching for an Ha nebula around RX J1605. 3+3249 
(e.g., van Kerkwijk & Kulkarni 2001a). A single definitive 
measurement would be the distance, which would deter- 
mine the areas of the X-ray and optical emission regions, 
but a parallax measurement with HST would require a 
significant investment of observing time. 
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Table 2 



Summary of Optically Detected Isolated Neutron Stars 



Source 


Period 


kT a 




log/A-//v c 


Optical 


References 




(s) 


(eV) 


(mag) 




Excess 




RX J1856.5-3754 




61 


25.8 


4.4 


6 


1,2,3,4 


RX J0720.4-3125 


8.39 


81 


26.8 


4.6 


6 


5,6 


RX J1308.6+2127 


10.31 


91 


28.7 


5.0 


5 


7,8,9,10 


RX J1605.3+3249 




92 


27.1 


4.4 


14 


11,12 



a Temperature of the best-fitting blackbody. 

b V-band Vega magnitude, either measured or interpolated. 

c Absorption-corrected bolometric X-ray-to-optical flux ratio, assuming that the X- 
ray emission is a blackbody. The V-band flux is computed according to fy = 
10 -(v+n.76)/2.5 ergs s -i cm -2 ; fo ii owing Bessell, Castelli, & Plez (1998). 

d The ratio of the observed F-band flux to the extrapolated X-ray blackbody flux at 
5500 A. 
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